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Алгоритм ранжирування атрибут!в для виявлення 
б1омаркертв у даних генно! експреси 


В статье рассматривается алгоритм ранжирования генов, полученных с использованием технологии 
микрочипов. Вектор рангов рассчитывается путем проведения классификаций случайных выборок из 
анализируемого набора данных. На каждой последующей итерации алгоритма ранг генов, участвующих в 
успешной классификации, повышается. В отличие от ранее используемых подходов, алгоритм позволяет 
повысить обобщающие свойства классификационных моделей за счет построения сбалансированных 
обучающих выборок, а также учесть информативность комбинации генов путем оценки их подмножеств. 

Ключевые слова: ранжирование признаков, биомаркер, классификация, генная экспрессия. 
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У статт! розглядаеться алгоритм ранжирування генв, отриманих з використанням технологй мкрочипв. 
Вектор ранмв розраховуеться шляхом проведення класиф!кашй випадкових виб!рок з анал1зованого 
набору даних. На кожн подальшй 1тераци алгоритму ранг генв, що беруть участь в успиинй класиф!кацй, 
шдвищуеться. На вдм1ну вд ранйпше використовуваних щдходв алгоритм дозволяе шдвищити уза- 
гальнювальн! властивост! класифлкащйних моделей за рахунок побудови збалансованих навчальних 
вибтрок, а також врахувати 1нформативн!сть комб1нацй генв шляхом ощнки 1х Шдмножин. 

Ключов! слова: ранжирування атрибут, б1омаркер, класифлкащя, експресля генйв. 


Введение 


Улучшение диагностики онкологических заболеваний является одной ИЗ перво- 
степенных задач, решаемых в области здравоохранения. Поставленный диагноз яв- 
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ляется основой для проведения и разработки курса терапии с целью достижения 
максимальной эффективности лечения с одновременной минимизацией токсичности. 
В последнее время большое внимание уделяется исследованиям генетической инфор- 
мации раковых клеток для уточнения существующих и выявления новых подтипов 
(классов) онкологических заболеваний на молекулярном уровне, что будет способ- 
ствовать эффективности последующего лечения. Технология микрочипов представляет 
собой новый класс биотехнологий, она позволяет измерять одновременно уровни 
экспрессии тысяч генов, что дает богатый материал для более полного понимания 
различий между подтипами опухолей на молекулярном уровне и, следовательно, для 
осуществления более достоверной их классификации. Данные, полученные с помощью 
технологии микрочипов, используются для определения генов с повышенной или 
пониженной экспрессией подтипов онкологических заболеваний, а также с целью 
выявления новых биомаркеров для клинической диагностики и прогноза [1], [2]. К со- 
жалению, извлечение зависимостей из такого рода данных затрудняется наличием 
несоответствия между количеством случаев (обычно менее ста) и количеством генов 
(обычно десятки тысяч), что приводит к переобучению, т.е. построению классифика- 
ционной модели с низкими обобщающими свойствами (низкая точность на незави- 
симой тестовой выборке) [3]. Кроме того, полученные данные, как правило, содержат 
как ошибки измерений, так и систематические ошибки, что может значительно повлиять 
на точность классификации. Поэтому классификации, так же как и кластеризации, 
должна предшествовать предобработка признакового пространства, а именно осуществ- 
ление отбора признаков (генов) путем ранжирования или сокращения их количества [4]. 

Методы отбора признаков для задачи диагностики онкологических заболеваний, 
как и в общем для всех задач классификации, можно разделить, в первую очередь, на 
две группы: фильтровочные и упаковочные методы [5]. Фильтровочные методы по- 
зволяют удалять неинформативные признаки согласно общим характеристикам 
данных. Они более просты в вычислительном плане и не привязаны к конкретному 
алгоритму классификации. Во встроенных методах выбранный классификатор участвует 
в построении подмножества признаков, тем самым, позволяя повысить точность 
классификации. Недостатком встроенных методов является высокая вычислительная 
сложность в связи с необходимостью на каждой итерации отбора применять классифи- 
кационный алгоритм. Для подмножества информативных генов, отбираемых из дан- 
ных генной экспрессии, применяются как фильтровочные методы, например, стандарт- 
ные параметрические тесты ({-тест [6]), непараметрические тесты (знаково-ранговый 
тест Вилкоксона [7]), так и упаковочные [8], [9]. Однако, зачастую, у разных авторов 
подмножества отбираемых генов различаются, что может быть результатом как не- 
соответствия размерностей анализируемых данных, так и игнорирования зависимос- 
тей между генами. Таким образом, на сегодняшний день не существует однозначных 
рекомендаций по выбору методов отбора подмножества информативных генов [10]. 

Цель данной статьи — разработка нового алгоритма отбора признаков, 
который позволяет преодолеть недостатки как фильтровочного, так и упаковочного 
подходов, и основан на ранжировании генов путем классификаций повторных 
случайных выборок из исходного набора данных. Предлагаемый алгоритм обладает 
следующими преимуществами: позволяет избежать переобучения за счет формирова- 
ния обучающих матриц меньшей размерности, с преобладанием количества случаев 
над количеством генов; позволяет учитывать информативность комбинации генов 
путем оценки подмножества признаков классификационным алгоритмом. 
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В статье представлены результаты тестирования предложенного алгоритма на 
наборе данных по лейкемии, оценена биологическая значимость подмножества из 
20-и генов с наибольшим рангом и проведен сравнительный анализ с результатами, 
полученными другими авторами. 


Описание предложенного алгоритма 


Пусть Хх мхм Представляет собой матрицу исходных данных экспрессии генов, 


где М - количество признаков (генов), а № -— количество случаев или объектов 
данных, вектор х, =(х,,х»,...,х”),1=ЪМ обозначает значения гена для № случаев 
или генный профиль. Каждый объект данных имеет метку класса, и в нашем иссле- 


довании соответствует одному из двух классов. Алгоритм автоматически масштаби- 
руется на обработку данных с большим количеством классов. Обозначим через 


№,1=1,2 количество объектов 7-го класса, тогда № = №, + М.. 

Алгоритм осуществляет ранжирование признаков путем классификации пов- 
торных случайных выборок из матрицы исходных данных экспрессии генов Ху, у. 
При этом выполняется процесс двойного отбора, т.е. случайным образом выбирается 
т<М генов и и<М объектов данных. Количество отбираемых объектов данных 
определяется параметром р, а количество генов такое, что т<и/2. Такая схема от- 


бора позволяет на А-Й итерации формировать сбалансированную матрицу У“ 


тхи > 

которую, согласно исследованиям авторов [3], менее вероятно переобучить. Коли- 

чество итераций К выбирается таким образом, чтобы каждый ген анализируемого 

набора данных отбирался определенное количество раз, достаточное для адек-ватной 

оценки его ранга. Дополнительным критерием останова итерационного процесса, 

кроме количества итераций, является значение коэффициента корреляции вектора 
рангов на текущей и предыдущей итерации. 

На каждом последующем шаге итерационного процесса, т.е. через определен- 

ное заданное количество итераций №, формируется выходная матрица Е. , с коли- 


чеством строк, равным исходному числу генов (1=1, М), и четырьмя столбцами. 
Значение первого столбца Т^ (1=1, М ) соответствует общему количеству раз, когда 


ген 7 включен в выборку для анализа после очередных А 


„ Итераций. Значение второго 
столбца 5" (1=1, М) соответствует общему количеству раз, когда ген { отобран в 


состав успешной модели классификации. Под успешной моделью подразумевается 
модель, имеющая ошибку классификации меньше предварительно заданного значения 


параметра &. Значение третьего столбца Р" =5*/Т^ (1=1,М ) соответствует предска- 
зательной способности 1-го гена и в четвертом столбце № (1=1,М) содержится 


значение ранга /-го гена, рассчитанное на основе Р^ после А„„, итераций. Гены с боль- 


Пег 
шим значением Р^ являются более информативными и, следовательно, имеют более 
высокий ранг. 

Каждая итерация алгоритма А =1, К состоит из выполнения следующих шагов: 


1. Для построения матрицы для анализа У, на К-й итерации применяется про- 


тхи 


цедура двойного отбора строк и столбцов матрицы из Х»„х. Матрица У, состоит 


тхи 
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из т случайным образом отобранных генов из всего исходного множества М, и, 
случаев из множества М, исходных случаев первого класса и и, случаев из множества 
№ исходных случаев второго класса, таких что и / М =и, / №, =В, п-+и, =п и 
т < п. Таким образом формируется обучающее множество, состоящее из и случаев, 
и тестовое множество, состоящее соответственно из №-—и случаев. 

2. К обучающей выборке, заданной матрицей У 


тхи > 


применяется классифика- 


ционный алгоритм. В качестве классификационного алгоритма используется метод 
сжимающихся центроидов [11]. Процедура перекрестной проверки используется для 
поиска оптимального классификатора, который имеет минимальную ошибку на 


обучающей выборке У*,. Использование процедуры позволяет определить пороговое 


тхи 
значение ДЛ”, соответствующее наименьшей ошибке классификатора, полученной с 
использованием наименьшего подмножества генов х,,х,,...,х,, где [< т. Полученная 
классификационная модель С“ = /(х,х,,....х,) верифицируется на тестовом множес- 
тве из М-и случаев. Ошибка классификации е рассчитывается следующим образом 
о ПРЕРМ 
ыы №М-п 
где ЕР — ошибка первого рода, Е№ — ошибка второго рода. 


В случае, если для всех возможных пороговых значений Л’ и подмножеств 
генов х,х,,....х,, [<т минимальная ошибка классификатора, построенного на обу- 


> 


чающей выборке, превосходит значение &, то верификация на тестовом множестве 
не выполняется и осуществляется переход к шагу 4. 

3. Если на тестовом множестве ошибка классификационной модели на тесто- 
вом множестве е^ < @, т.е. меньше, чем предварительно заданный порог, то модель 
принимается как успешная и выполняется модификация выходной матрицы Ё\,. 


Для каждого х, гена анализируемой на К-й итерации матрицы У . определяется 


тхи 


соответствующий ген матрицы Е„„. и значения столбцов 1^, 5 и Р" рассчиты- 
ваются следующим образом: 


2°т 


в: о Е 5+1, хе (м,х,....%) Ре 114 


1 


к-1 КЬ- 
А о, ен 


т 


4. Если на тестовом множестве ошибка классификационной модели г“ >а, 
классификационная модель считается неуспешной И отобранные на К-й итерации 
гены являются неинформативными. Модель в этом случае можно определить как 


переобученную на обучающей выборке. Значения столбцов Т^, 5’ и Р" выходной 
матрицы Е„„. для каждого гена из отобранного подмножества (х,,х,,...,х„) рас- 
считываются как 


к-1 
Т + хЕ(х,х % 
И 1 2 1 а. К _ © К _ © К 
=] те в =, У Й 
о нак 


Значения рангов определяются путем сортировки прогностических значений ге- 


т, ыы 
нов Р', 1=1М в порядке убывания и заносятся в столбец А, 1=1,М матрицы Ед. 
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Каждые [е итераций оцениваются изменение в ранжировании генов путем рас- 


тег 
чета коэффициента ранговой корреляции Спирмена между векторами К и №` на 
К-йи К* -Й итерациях, где А*=А-К,,: 


Пег 
Е ЕО х 
К _ 7 2 7 12722. **? 


А КЬ-1 К К К 
ри о о ФР 
р Ох 


т 


Значения рангов определяются путем сортировки прогностических значений 
генов Ру 1=1, М в порядке убывания и заносятся в столбец К ‚ ТБМ матрицы Ед. 


Каждые [е итераций оцениваются изменения в ранжировании генов путем рас- 


тег 
чета коэффициента корреляции Спирмена между векторами рангов А и А“ на К-й 
и К* -й итерациях, где К*=А-К„: 


Пег 


Алгоритм прекращает свою работу в случае, если коэффициент корреляции 
г>у (например у = 0.99 ), т.е. когда ранговая последовательность стабилизируется. 


Р езультаты экспериментов 


Для тестирования предложенного алгоритма ранжирования признаков и про- 
ведения сравнительного анализа результатов в нашем исследовании использовался 
набор данных пациентов с лейкемией [12]. Набор данных включает два типа острой 
лейкемии: острая лимфобластная форма (АГЛ.) и острая миелоидная лейкемия (АМГ.. 
Набор состоит из 47 случаев лейкемии АПЛ, (38 случаев В-се| лейкемии и 9 случаев 
Т-сей лейкемии) и 25 случаев АМГ, которые характеризуются экспрессией 7129 генов. 
Согласно протоколу, описанному в [13], были выполнены следующие шаги предоб- 
работки данных: пороговая обработка с наименьшим значением экспрессии, равным 
100 и наибольшим — 1600; фильтрация с исключением генов со значениями экс- 


прессии, удовлетворяющими одному из следующих условий: и < 5 или 


(плах— 1011) <500, где тах — максимальное значение и шш — минимальное значение 


экспрессии гена; логарифмическая трансформация значений экспрессии генов. Пос- 
ле выполнения предобработки было отобрано 3571 генов. Набор данных можно 
скачать по ссылке [14]. 

Для проведения эксперимента были выбраны следующие параметры: коли- 
чество итераций К =300000, пороговое значение ошибки классификации &@ = 0.2, 
значение доли случаев, случайным образом отбираемых из исходного набора на 
каждой итерации Д=0.7, количество итераций между выполнением расчета кри- 


терия останова К 


Пег 


=20000 путем оценки коэффициента корреляции с пороговым 


значением г = 0.09. Каждые А‚„ оценивалась выходная матрица ЁЕ,„.. В табл. 1 пред- 


[ег 
ставлены 20 генов с наибольшим рангом, полученные в результате эксперимента. 
Для каждого из 20 подмножеств, начиная с одноэлементного, состоящего из 
гена с наивысшим рангом, и далее путем добавления последующего гена в ран- 
жированном списке, была построена классификационная модель и оценена ошибка 
классификации всего исходного набора данных. В качестве наиболее информа- 
тивного подмножества, которое может рассматриваться как набор потенциальных 


биомаркеров заболевания, было отобрано подмножество с наименьшей ошибкой 
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классификации. На рис. | представлены значения ошибки классификации, соответ- 
ствующие количеству отобранных генов от 2 до 20. Согласно полученным значениям 
ошибки, рассчитанным с использованием процедуры перекрестной проверки на всем 
исходном наборе данных, в качестве биомаркеров заболевания можно выбрать 
подмножество, состоящее из четырех генов (С433, Хухш, АРГ.Р2, С$ТЗ), с точностью 
классификации 98,6%. Необходимо отметить, что подмножество из 10 генов дает 
несколько меныпую ошибку классификации, однако требует определения значений 
экспрессии большего количества генов. 


Таблица 1 — Описание генов с наибольшим рангом 


[№ Обозначени И № Обозначение Е 
е гена гена 
| | Са33 Рекомбинантный | рмтТ Терминальная диоксинуклеотид- 
белок человека трансфераза 
2 | Хухш Циксин 12 | МАВСК$Г1 | МАВСК$5-подобный протеин 
8 | АРГР2 Предшественник 13 | СР79А Полипептид 1э-альфа (иммуно- 
бета-амилоида глобулин ассоциированная альфа) 
Ц | С5ТЗ Цистатин С 14 | СТЗА Катепсин А 
5 | МОУ и 15 | СЗТА Цистатин А 
трансферазы 
6 | СТО Катепсин О 16 | ЗЕВРИМВ1 Серпин В1 
7 | СЕ Адипсин 17 Г ЕАН Фумарилацетоацетат гидролаза 
Пропердин — глобулярный белок 
8 | ССМОЗ Циклин 03 18 | СЕР сыворотки крови человека и 
млекопитающих животных. 
Мембранный о цепи иммуноглобулинов. 
ВпОК, о —. который 
9 | СО63 гликопротеин из | 19 | УРВЕВ1 ны не 
ее стволовая клетка), но который 
семейства 
активируется во время ранних 
тетраспанинов в 
стадий развития В-[утрросуе. 
10| ТСЕЗ а 20 | ЗРТАМ! Спектрин альфа-цепь 
транскрипции 
21 | МРО Миелопероксидаза 


Ошибка классификации (количество ошибок) 


0.0 0,5 1.0 1,5 20 25 30 35 40 45 50 55 60 
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Рисунок 1 — Ошибка классификации для количества генов от 2 до 20 
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Сравнительный анализ результатов, полученных с использованием предложен- 
ного алгоритма, показал его эффективность при определении наиболее информа- 
тивных генов, Т.к. большее их количество присутствует в списках информативных 
генов других авторов [15-18], а сравнимая точность классификации (с использо- 
ванием процедуры перекрестной проверки) достигается с использованием меньшего 
их количества (табл. 2). 


Таблица 2 — Результаты сравнительного анализа различных методов отбора 


признаков 
Точность классификации | Количество генов 
НВЕ метод [15] 97.2% 4 гена 
Метод [16] 98.6% 132 гена 
Метод [17] 98.6% 3 кластера генов (размер одного 
кластера от 1 до 23 генов) 
Предложенный алгоритм | 98.6% 4 гена 


Анализ биологической значимости отобранных генов 


Для того чтобы оценить биологическую значимость полученных результатов 
ранжированное подмножество из 20 наиболее информативных генов было проана- 
лизировано с использованием инструмента анализа биологических процессов шпаерВ 
(Бер://\у\у\ .шпае4Ь.са). Была выявлена их потенциальная функциональная связь с 
фенотипом заболевания и биологические пути, в которых они участвуют. Значимость 
участия комбинации генов в определенном биологическом процессе оценивалась путем 
расчета соответствующего р-уровня. Согласно полученным результатам анализа два 
(адипсин СЕР и комплементарный фактор пропердин СЕР) из 20 генов участвуют в 
альтернативном пути системы комплемента, подавлении инфекционных агентов (р- 
уровень = 0.00006). В работе [18] адипсин был также выбран в качестве биомаркера 
подтипа лейкемии, и была описана его роль в дифференциации миелоидных клеток. 
Два гена (СОЗЗ и ОМТТ) из 20 связаны с процессом дифференцировки 
гемопоэтических клеток (р-уровень = 0.00429). СОЗЗ кодирует мембранный протеин 
клетки и сильно экспрессирован на поверхности лейкемических бластов, приме- 
няется для дифференциации миелоидных и лимфоидных клеток [19]. ОМТТ 
(деоксинуклеотидтрансфераза) экспрессируется в популяциях нормальных и злока- 
чественных пре-В и пре-Т лимфоцитах во время ранних стадий дифференцировки. 

Два гена (циклин ОЗ и СО79А) участвуют в биологическом процессе клеточ- 
ного цикла и сигнального пути В-клеточого рецептора (р-уровень = 0.01489). Также 
два гена циклин ОЗ и циксин участвуют в процессе фокусной адгезии (р-уровень = 
=0.0225). В работе [20] циксин был отобран как один из биомаркеров фенотипов 
лейкемии. Уровень его экспрессии играет важную роль в дифференцировке подтипов 
лейкемии, однако его прямая связь в гемопоэзе не выявлена. Согласно последним 
исследованиям циксин-кодируемые протеины могут регулировать транскрипцию гена 
путем взаимодействия с фактором транскрипции. Таким образом, наиболее инфор- 
мативные 20 генов с наибольшим рангом, полученные с использованием предложенного 
алгоритма, функционально значимы как для процесса дифференцировки гемопоэти- 
ческих клеток, так и патогенеза онкологических заболеваний. Многие из отобранных 
нами генов отобраны также в качестве биомаркеров для распознавания подтипа 
лейкемии в ряде других исследований [18-20]. 
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Оценка сходимости алгоритма 

Для определения оптимальных прогностических значений каждого гена необ- 
ходимо оценить все возможные выборки т случаев из исходной матрицы данных 
Хуи»х, а именно Су, = (М = 1 (т = 1) (М = т)! различных комбинаций, что является 
сложным в вычислительном плане с количеством вычислений, пропорциональным 
О(М!). Используемый нами подход является приближенным, эвристическим и по- 
лучает оценки прогностических значений генов путем повторного анализа случайных 
выборок 77 случаев из исходной матрицы данных Х»„„у. Следовательно, необходимо 


проанализировать сходимость предложенного алгоритма и способность получать 
достоверные оценки ранжированных прогностических значений. Для этого мы про- 
анализируем вектор рангов для различного количества итераций алгоритма. Как 
показано на рис.2 для 20 генов, расположенных вверху ранжированного списка, 
значения рангов на 20 000 и 40 000 итерациях достаточно сильно отличаются. 


—$Ф— М23197_а+ 
—\—Х95735_а4 
—__ [09209_$_а+ 
—>— М27891_а4 
—Ж— 146499 а 
—®— \63138_а+ 
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Рисунок 2 — Изменение рангов 20 генов в зависимости от количества итераций 


Однако, после приблизительно 200 000 итераций значения рангов стабилизи- 
руются. Таким образом, предложенный алгоритм сходится к близкому к оптимальному 
решению, и может достаточно эффективно осуществлять ранжирование признаков, 
используя сравнительно небольшое количество итераций. 


Выводы 


Предложенный алгоритм ранжирования признаков позволяет выделять наи- 
более информативные гены в данных генной экспрессии, причем стабильность значе- 
ний рангов отдельных генов обусловлена многократной оценкой выборок из исходной 
матрицы данных, что помогает избежать переобучения, и способствует получению 
несмещенных оценок вектора рангов. Согласно алгоритму на каждой итерации выпол- 
няется классификация случайным образом сформированной обучающей выборки с 
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верификацией классификационной модели на тестовой выборке. Точность классифи- 
кации является показателем прогностической способности отобранных на очередной 
итерации комбинации генов. Выходная матрица количественно фиксирует как прог- 
ностическую способность, так и ранг отдельных генов, которые модифицируются на 
каждой последующей итерации, приближаясь к оптимальному ранжированию. Крите- 
рием оптимальности вектора рангов является стабилизация его значений для отдель- 
ных генов, что оценивается с использованием критерия ранговой корреляции Спирмена. 

Предложенный алгоритм протестирован на наборе данных по лейкемии и про- 
ведена оценка сходимости алгоритма на примере 20 генов, расположенных вверху 
ранжированного списка. Оценена биологическая значимость исследуемого подмно- 
жества отобранных генов, которая позволяет сделать вывод об их тесной связи с 
процессами, происходящими в лейкемических клетках. Таким образом, отобранные 
гены с наивысшим рангом не являются результатом случайного отбора. Сравни- 
тельный анализ с работами других авторов по исследуемому набору данных показал 
преимущество предложенного нами алгоритма, а именно в качестве результата ото- 
брано наименьшее подмножество из четырех биомаркеров, обеспечивающих такую 
же или лучшую точность классификации. 
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КЕ5ОМЕ 
№. № №оуо5е[ота, 1. Тот 


Аогийт о} Кешиге Капюиз ог Втотагкег О15соуегу т Сепе 


Ехргезюп Ба 

п Фе рарег Фе гапКто а1еогИт о эепе ехргез$1оп даа оМаштей Нот писгоаггау 
теазигетеп $ 15 сопз14егеа. ТВе ргорозе4 а!еогит аПоууз тю зейесе Фе поз шогтануе 
сепез Бу гапКшо, \УПеге Фе за Шу оЁ Фе гапк$ ог шагуиа| эепез 15 епзиге4 Бу 
езитайоп оЁ тире затр!ез Нот шШа! даёа таблх. ЗисВ ап арргоасВ Бе]рз ю ауо!4 
оуег по ап епаез$ Ше ипЫазе4 езита{е оЁ Ще уесг оф гапКз. А| еасп Цегайоп Фе 
с1а51ИсаНоп то4е| сопзиасе4 оп Фе гапдоп у эепегае4 гаште затр/е 1$ уемйе4 оп 
{Бе {е5{ затр!е. ТБе с1аз5ШсаНоп ассигасу 1$ 1е шсаюг оЁ е ргоспозйс аБПиу оЁ Фе 
шагу ча! оепез. АЙег фе зиссез$ В с1азШИсайоп Ше гапК$ оГ Фе рагистрайпе сепез 
Бесоте №1орег, теапуИе е зеагсВ г Фе орйта| гапК1е 15 регогтеЯ по{ Рог еасв 
шагу а! эепе, Ба Фе у’вое сот Ыпайоп. ТНе опёриё таблх 15 шоШе айег рге- 
да&егтте4 питЪег оЁ ЦегаНопз, геол его Бо Ще ргоспо$Нс аБШиу ап4 е гапК$ ог пе 
шагу1Аа| оепез. ТВе за Шу оЁ Фе гапК уесюг зегуез аз ап оритаШу стйепоп апа 1$ 
езитае4 Бу сотрийпе Зреагтап сотге]аНоп соеРЯслепе Бебмуееп Ше сигтеп ап4 ргеу1юиз 
гапК ог4ет. 

ТВе ргорозе4 а!еогИфт Ваз Бееп 1е${е4 оп фе 1еиКепта даёазе{ ап4 Из сопуегоепсе 
13 апа|у7е4, сопуАегте Фе \епбу 'юр-гапКе4 сепез. Те апа[уз15 оЁР е ЫМо]оз1са| 
юпИсапсе ог фе шуезИозме оепе зибзеЁ аПоу’з ю сопйгт И5 об\у1юи$ КпсНопа| 
ге]еуапсе 1ю Ше рБепобуре И рге1с{$ ап4 Ше ргосеззе$, факте расе ш Ве |епКепис се. 
[ аззигез ШаЁ Фе юр-гапКе оепез аге Ше Му ипПКеу тю Ъе з@есцеа Бу сВапсе. ТВе 
сотрагайуе апа[у$15 оЁ Фе 4еуеоре а]еогИйт оп Фе |епКепиа афазеё зВо\уз И5 
адуащасе оуег апа|о2$, побау фе зейесе@ зеё ог МотагКегз 1$ зтаПет, соп$15Я по оЁ 
Гопг оепез, УсВ ргоу14е зипПаг ог ШоВег с1аз$1Исайоп ассигасу ап@ ргезегуе Фе №121 
с1аз1ИсайНоп ассигасу. 
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